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Wireless objects equipped with multiple antennas are able to simultaneously transmit or receive 
multiple packets by exploiting the channel's spatial dimensions. In this paper, we study the benefits of 
such Multiple Packet Transmission/Reception (MPTR) approach, when it is used in combination with a 
CSMA/CA MAC protocol for fully interconnected networks, addressing the interactions between the two 
mechanisms and showing the performance gains that can be achieved. To this end, a new MAC protocol, 
called MACMA (MAC protocol for Multiple Antenna wireless objects), is introduced and formulated as a 
very basic MAC protocol that captures the fundamental properties of a CSMA/CA channel access protocol 
supporting MPTR. Using the MACMA protocol as a reference, a new mathematical model for CSMA/CA 
MAC protocols supporting MPTR is presented for the case of non-saturated traffic sources with finite 
buffer space. Simulation results show that the analytical model is able to accurately characterize MACMA's 
steady-state behavior for different number of antennas and different traffic loads, providing a useful tool to 
understand the performance gains that the proposed scheme can achieve. 
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Modeling a CSMA/CA MAC Protocol with 
Multi-packet Transmission /Reception in 
Non-saturated Conditions 

1 Introduction 

MIMO techniques allow for a more efficient use of transmission resources by making 
the channel's spatial dimension available, using simple and effective signal processing 
techniques II). This new space dimension can be used to enhance coverage and reliability (spatial 
diversity) or to boost the system throughput (spatial multiplexing) Q. In this paper, we will be 
focused on the spatial multiplexing feature, which can be used either to transmit or to receive 
multiple packets at the same time, referred to as Multi-Packet Transmission (MPT) and Multi- 
Packet Reception (MPR), respectively When both options are available at the same node, we 
say that it has MPTR capabilities. MPT (MPR) can be seen as a packet-based extension of the 
Multi-user MIMO downlink (uplink) concept, which is also known as Space Division Multiple 
Access (SDMA). 

Although the benefits and drawbacks of MPTR in point to point links are well-known, there is 
still a lack of results focusing on how these MPTR capabilities can be efficiently combined with 
random access MAC protocols, in order to optimally exploit the performance gains that this 
combination may offer. On the one hand, using MPTR will reduce the number of transmission 
attempts and therefore decrease the collision probability in a contention scenario. On the other 
hand, in order to make the MPTR work properly, the MAC protocol needs to be modified, e.g., 
to include the necessary procedures for estimating the channel state between the transmitter 
and the receiver, required for isolating the multiple spatial streams, or to feed the CSI (Channel 
State Information) back to the transmitter in those schemes that need it. There will also be 
a need for extra acknowledgments to confirm the reception of all spatial streams. Both of 
these modifications will impose extra temporal overheads, the negative impact of which on the 
system performance needs to be quantified and contrasted with the performance boost gained 
by enabling the transmission of multiple streams. 

This paper intends to provide insight on the interactions between CSMA/CA random access 
MAC protocols — taking the IEEE 802.11 Distributed Coordination Function (DCF) [3] as a basis — 
and the MPTR scheme. We study and characterize these interactions by introducing a new 
protocol that combines these two mechanisms, and evaluating the performance of the overall 
system using an analytical model, as well as simulations. The presented mathematical model 
provides a formal representation of the new MAC protocol, helping to understand its behavior, 
and is a suitable tool for finding robust parameter values that yield to a close-to-optimum 
performance. 

In this context, the main contributions of this paper are: (1) The introduction of a new MAC 
protocol, called MACMA (MAC protocol for Multiple Antenna wireless objects) that comprises 
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the fundamental operations required for an MPTR CSMA/CA-based random access protocol, 
including those required for channel estimation and a new ARQ (Automatic Repeat reQuest) 
protocol for retransmission of the packets that are lost due to errors or collisions; (2) An analytical 
model for the MACMA protocol, based on combining a batch-service queueing model with a 
CSMA/CA MAC protocol model in non-saturation conditions, showing that this combination 
provides a suitable path to model packet-based multiple antenna access protocols; and (3) The 
performance evaluation of the MACMA protocol, providing insights on how its parameters have 
to be configured and showing the performance benefits of such kind of protocols, as well as 
confirming the accuracy of the presented model. 

The paper is structured as follows. Section [2] summarizes the current state of the art on MAC 
protocols and analytical models for multi-antenna nodes. Section [3] presents the considered refer- 
ence scenario, i.e., a single-hop network where N nodes compete for the channel. Section [4] details 
the MACMA protocol. In Section [5} the analytical model is presented. The analytical model 
validation and performance results are discussed in Section |6j Finally, the main conclusions of 
the paper are summarized in Section 

2 Related Work 

CSMA-based MAC protocols with MPT capability have been considered for both Ad-Hoc net- 
works BU, 13 and WLANs 0, due to their ability to improve the network performance, 
while keeping both the simplicity and the efficiency that random access MAC protocols possess. 
In the aforementioned works, the MAC protocol is extended to both coordinate the multiple- 
packet transmissions and to provide the required CSI to the transmitter, the potential receivers, 
or both. Even though these proposals have been evaluated through simulations, only in a few 
cases analytical models have been presented @, 0. Moreover, the queueing dynamics in such 
protocols have not yet been analytically evaluated, which is one of the key contributions of our 
paper. 

In HI, a MAC protocol for WLANs supporting MPT at the Access Point (AP) is presented. 
The RTS/CTS handshake procedure is extended to both coordinate the transmission of multiple 
packets from the AP to different STAs, and to provide the AP with the required CSI. An analytical 
model of the presented protocol in non-saturation conditions is introduced, although queuing is 
not considered. The non-saturation case is approximated through normalizing the transmission 
and collision probabilities by the probability that the AP has packets present in its queue, which 
is similar to the approach we use in this paper. However, they assume that once the AP's queue is 
not empty, it has as many packets as the number of antennas present, which may not necessarily 
be the case. 

In 0, the performance of MPT in the downlink for the upcoming IEEE 802.11ac standard 
Il8l is evaluated. An accurate physical-layer and channel model is provided when a zero-forcing 
precoding scheme [93 is used to create the multiple spatial streams. To evaluate the system 
performance, an analytical model in saturation conditions is presented, where it is assumed that 
the AP can always transmit as many packets as the number of antennas it has. 
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In our work, in order to analytically model the MPT capability in non-saturation conditions 
and focus on the queueing dynamics, we consider the use of batch-service queueing models lHUl , 
which are currently receiving a lot of attention due to their application in modeling various 
features of communication systems, such as packet aggregation schemes on a temporal scale 
|TT|. Examples of use in specific technologies are |12|, 11131 , |T4| focusing on the performance 
analysis of temporal aggregation in Wireless Local Area Networks (WLANs). 

It should be noted here that, traditionally, CSMA/CA-based MAC protocols with MPR capa- 
bility have been mainly considered for WLANs, where an AP equipped with multiple antennas 
is able to receive multiple packets simultaneously Ifl5ll , Ifl6l , usually transmitted from different 
nodes. In MACMA, all nodes are also equipped with multiple antennas and have MPR capability, 
as they have to receive and decode all the transmitted packets from a transmitter equipped with 
multiple antennas. However, in MACMA simultaneous transmissions from multiple transmitters 
are treated as collision, and MPR capability is only used for receiving multiple packets from a 
single transmitter. 



3 System Model and Assumptions 

A single-hop network with N nodes is considered. Each node is assumed to be within the 
transmission range of the others and close enough to have negligible propagation delay. In 
addition, regardless of the specific positions of the transmitting and receiving nodes, it is as- 
sumed that for every transmitted packet, all the receivers observe the same average SNR. This 
assumption allows us to consider all nodes in the same collision domain, which results in all 
transmitted packets having the same Packet Error Probability (PER ). This greatly simplifies 
the model formulation, without affecting the MPT queuing model, which is one of the main 
contributions of the paper. 



3.1 Node Operation and Link Layer 

Each node is equipped with M antennas and a finite-buffer of K packets, as shown in Figure 
[TJ The MPTR capability by using the MIMO Spatial Multiplexing mode is considered, allowing 
the transmission of multiple packets simultaneously from a single sender to a single or up to M 



4 



destination nodes. A similar approach is considered in [T7|, where the sending node allocates 
the transmit antennas independently and based only on the destinations of the selected packets 
and, therefore, a single user can receive none, one or more packets at each transmission. Packets 
arrive to each node following a Poisson process with rate A. Each packet has a fixed length of 
L bits and can be directed to any of the other nodes. 

At each transmission, between s min and s max packets can be transmitted simultaneously through 
the s created spatial streams. Both s min and s max parameters can take values between 1 and M, 
with s min < s max . Packets are selected for transmission following a First-In First-Out (FIFO) 
policy, regardless of their destination(s). Packets depart from the queue in batches, called space- 
batches, of size s(q) = max{s min , min{g, s max }} packets, where q G [0, K\ is the queue occupancy 
just before the new transmission is scheduled. If q < s min , the link layer remains inactive until 
enough packets have arrived. Choosing a high value for s m i n reduces the number of transmission 
attempts on the channel, thus reducing the probability that a transmission results in a collision 
but at the cost of a larger average waiting delay, specially at low traffic rates. The s max value has 
to be adjusted considering the channel state. Under good channel conditions (high Signal-to- 
Noise Ratio (SNR)), s max can be increased towards M as the MACMA can benefit from a larger 
number of parallel transmissions with fewer transmission errors. However, at low SNR, high 
Smax values will result in a high packet error rate (PER ), hence limiting the system throughput. 

The channel access is governed by the MAC protocol presented in Section HI It is called 
MACMA (MAC protocol for Multiple Antenna wireless objects) and is designed to exploit the 
presence of multiple antennas in the wireless objects to increase the network performance in 
terms of throughput, delay, and number of nodes that can be supported in the network. 

3.2 Operation of Other Layers 

Since the focus of this paper is on the link-layer (transmission queue and MAC protocol) and 
in particular on an analytical model for computing the impact of link-layer parameters on the 
system performance, an idealized channel, physical layer, and source traffic model are assumed 
in the system model. 

3.2. 1 Channel Model 

Multi-packet transmission allows for transmitting multiple packets simultaneously from the 
same node to one or more destinations. 

A quasi-stationary channel with flat fading coefficients changing independently from trans- 
mission to transmission (MIMO Rayleigh channel) is considered. 

3.2.2 Physical Layer 

In transmission, the physical layer, after receiving a group of m selected packets from the link- 
layer, is responsible to (i) build the space-batch frame by mapping each scheduled packet to an 
antenna and {%%) encode and modulate the signal for its transmission over the channel. Equal 
transmission power of P t /m is allocated to each packet, where P t is the total transmission power. 
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Two transmission rates are considered; one for the physical layer headers, r p h y/ and the other 
for the MAC protocol headers and data, r. 

In reception, the physical layer performs the (i) transmitter /receiver synchronization, (ii) 
channel estimation, and (Hi) detection of the different spatial streams using a Zero-Forcing (ZF) 
detector 10, which are successively de-modulated and de-coded. The average received SNR at 
each node is £ and £o/m is the average received SNR at any receive antenna contributed by a 
single transmit antenna. 

Considering previous assumptions and a ZF detector, the post-processing SNR at each spatial 
stream is usually assumed independent of the other branches and follows a x 2 distribution with 
1 = 2- (M — m + 1) degrees of freedom, where M is the number of receiving antennas and m 
the number of simultaneously transmitted spatial streams [9J. Based on this, and assuming that 
a packet is received correctly only if the post-processing SNR is higher than a certain reference 
SNR, £ re f, the Packet Error Rate (PER) when the transmission contains m spatial streams is 
equivalent to the probability that the post-processing SNR, £, is smaller than £ ref , i.e., 




where is the cumulative distribution of the post-processing SNR. Note that, given a fixed 
number of antennas, the PER increases with the number of spatial streams transmitted, thus 
showing the existence of a reliability-capacity tradeoff [2j. Capture effect is not considered and 
in case of collision, all transmitted packets are assumed to be corrupted and cannot be decoded 
correctly. 

4 MACMA: A MAC Protocol for Multiple Antenna Wireless Objects 

The MACMA extends the IEEE 802.11 Distributed Coordination Function (DCF) |3l and, unless 
the contrary is stated, it shares with the DCF the channel access and backoff techniques, and 
the meaning and functionality of parameters such as the DIFS (Distributed Inter-frame Space) 
and SIFS (Short Inter-frame Space) parameters. 

4.1 Protocol Description 

At any given time, a node operating under the MACMA protocol can be in either of the following 
operation modes: idle, transmission, and reception. 

4. 1. 1 Idle Mode 

A node is in idle mode if it is neither transmitting nor receiving. From the idle mode, the 
MACMA will move to the reception mode if the channel is detected busy, or to the transmission 
mode if the number of packets stored in its queue reaches s min due to new arrivals. 
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4.1.2 Transmission Mode 

As soon as a node enters transmission mode, it builds a new space-batch using the first s(q) — 
max{s min , min{g, s max }} packets stored in its queue. It will then wait until the channel is detected 
free for a period equal to a DIFS after which it initiates a random backoff to avoid collision with 
transmissions from other nodes. When the backoff ends, the space-batch is transmitted and the 
transmitting node waits for the corresponding ACKs. 

When a transmission is scheduled, the backoff counter is randomly set to a uniformly dis- 
tributed value picked from the range [0, CW — 1], with CW the length of the contention window. 
The backoff counter is decreased by one unit at the beginning of every slot until it reaches 
(selecting a random value equal to means immediate transmission) and the space-batch is 
transmitted. During backoff, if the channel is detected busy, the countdown is frozen until the 
channel becomes idle§ again. 

A transmission can have three different outcomes: it can be successful, end in collision, or 
arrive collision-free but contain erroneous packets. We say that a packet is transmitted success- 
fully if its reception is acknowledged by the receiver, which means that it has neither suffered 
a collision nor contains errors. 

Figure |2] shows an example of the operation of MACMA for both cases of a collision-free 
transmission as well as a collision. As shown in the figure, a separate ACK is transmitted 
for every successfully received packet. The packets corresponding to a space-batch will be 
purged from the queue only when the reception of all packets in that space-batch have been 
acknowledged. When a packet is not correctly received, the same retransmission mechanism 
is applied regardless of whether the loss has been due to collision or transmission errors. 
In this case, MACMA applies a space-batch-oriented retransmission mechanism that at each 
transmission, sends only the remaining (unacknowledged) packets of the previous space-batch 
until all the packets are correctly received. This means that the newly admitted packets to the 
queue will not be transmitted within the retransmitted space-batches even if the number of 
retransmitted packets included in the space-batch is lower than s max . This can be observed in 
the fourth transmission of the second node in Figure |2l where only the unacknowledged packet 
from the previous transmission is sent. The motivation of this approach is that, the fewer the 
number of packets transmitted in a space-batch, the higher the chances that they do not suffer 
transmission errors and therefore the expected retransmissions due to errors is kept low. 

The time during which the channel remains busy due to a collision-free transmission or 
retransmission of a space-batch of size m, T c f (m), regardless of whether or not it contains errors, 
is given by 

Trf (m) = Tdata(m) + ACK time out (2) 

where Td a t a (m) is the duration of the data portion of the transmitted (or retransmitted) frame 
that includes m packets, and ACK timeout is the time spent waiting for and receiving the ACK 
for each successfully transmitted packet. Note that, to avoid confusion, we use the notation s 

1. The channel is detected idle only if it has remained idle during at least DIFS units of time. 
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Fig. 2: MACMA protocol: Successful Transmission and Collision Duration for s m j n = 1, s max = 2 and M = 2. 



to refer to the original space-batch size, and m < s, to refer in general to the number of packets 
involved in a transmission or a retransmission. 

It is important to note here that ACK timeout has to be set to its maximum possible value, i.e., 

ACK timeout = M ■ (SIFS + T ack ) (3) 

where SIFS is the Short Inter-Frame Space as defined in (3J and T ac k is the duration of an ACK 
frame. This is because, there can always be nodes that are not able to decode the space-batch 
due to errors and therefore, cannot know the number of packets that it includes. Furthermore, 
a node that does receive the space-batch successfully has no way of knowing whether all other 
nodes have been able to do the same. Since the MACMA requires that all nodes restart the 
channel contention at the same time after a transmission, all nodes will have to wait the total 
time required for receiving M acknowledgements. 

Similarly, in case of a collision, both the transmitter and the nodes that have overheard the 
transmission have to wait for a full ACK t i me0 ut in order to receive all possible ACKs. Additionally, 
as space-batches of different durations, carrying a different number of packets, can collide, the 
ACK timeout has to be initiated only after the transmission of the data part of the longest space- 
batch involved in the collision. This guarantees that all nodes will restart to contend for the 
channel at the same time. Therefore, the collision duration is given by 

T c (m) = T dat a(m max ) + ACK tim eout (4) 

where m max represents the longest space-batch involved in the collision. 

Note that since the transmitter does not have access to the CSI (Channel State Information), it 
is not able to decode the ACKs that are simultaneously sent from all the receivers. Nevertheless, 
if this was possible, the ACK timeout could be significantly shorter, resulting in an improvement in 
the overall system performance. However, this would be at the cost of having extra overheads 
necessary to feed the transmitter with the CSI from each receiver, e.g., by using a modified 
RTS/CTS mechanism. 

4.1.3 Reception Mode 

A node in reception mode checks the receiver address of each correctly received packet from the 
arriving space-batch frame. For each packet that is directed to it, it schedules the transmission of 
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Fig. 3: Space-batch structure as considered by the MACMA protocol 

an ACK following the order indicated in the Stream Control (SC) field of each data packet. This 
means that the i-th ACK is transmitted SIFS+(i— 1) (T ac k + SIFS) seconds after the reception of the 
space-batch. After the ACK timeout , all nodes move to the idle mode or remain in the transmission 
mode, depending on their queue status. 



4.2 Space-Batch and ACK Frames 

The basic space-batch and ACK frames defined by the MACMA protocol are depicted in Figure 
|3l They include the specific PHY and MAC headers required to support the described physical 
and link-layer functionalities. A detailed description of the headers is not needed for the present 
paper. Nevertheless, for the sake of completeness, we briefly outline their contents: (i) the Signal 
(SIG) field includes the number of training sequences that are sent and the antenna to which 
they correspond, (ii) the Frame Control (FC) includes the specific information about the entire 
transmission (i.e., type of frame (data or ACK), number of frames transmitted, etc.) and (Hi) the 
Stream Control (SC) includes the specific information for each frame inside the space-batch (i.e., 
sequence number, retransmission flag, etc.). 

The duration of the space-batch frame is variable and depends on the number of training 
sequences included. Note that for the specific case of m = 1, the training sequence is not required 
for the receiver as there is no need to detect multiple transmitted streams. Therefore, the duration 
of a space-batch frame is given by: 



Tdata(m) 



rphy r 
L p hy+L tt (m) 

r p hy 
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(5) 



where L phy , L tr (m), and L mac are the PHY header, the variable number of training sequences, 
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and the MAC header length, respectively, as shown in Figure [3] The ACK has a fixed duration 
equal to T ack = ^ + ^ seconds. 

5 Modeling the MACMA 

The MACMA model consists of two main parts: (1) A Batch-service queueing model used to 
capture the multi-packet transmissions, and (2) The CSMA/CA non-saturated MAC model used 
to capture the interactions among nodes when accessing the shared channel. The MAC model 
includes also the specific space-batch-oriented ARQ protocol. 

The MACMA model is solved numerically by applying a fixed-point approach. Given an 
initial set of parameters, the fixed-point algorithm iterates until the convergence condition is 
reached (i.e., all the parameters of the model stabilize to a specific value). Note that both the 
batch-service queueing model and the MAC model are coupled by cross-related variables such 
as the service time (X) and the probability to have packets ready for transmission (p). 

5.1 MPT Queueing Model 

An M/G^/l/K queue is used to model each node. In this notation, [s] = [s min ,s max ] indicates 
an interval, with s min and s max representing the minimum and maximum number of spatial 
streams that can be used at each transmission, respectively. Poisson arrivals of rate A and a 
general service time distribution are considered. The buffer has a size of K packets and no extra 
space is considered for the packets in service, and therefore, the packets included in a space-batch 
transmission remain stored in the queue until they are correctly received and acknowledged or 
discarded after the maximum retry limit is exceeded. 

As mentioned before, packets are served in batches of size s(q) = max{s min , min{g, s max }}, i.e., 

{^min Q ^ ^min 

q s min <q< s max (6) 
■Smax ^max — Q 

where q is the number of packets present in the queue after a departure. The scheduling of a 
space-batch takes place after each departure, as soon as enough packets become available in the 
queue. Here a departure refers to the moment at which all packets from the previous space-batch 
are acknowledged and purged from the queue. If just after a departure there are more than 
■Smin packets in the queue, the next space-batch will be scheduled immediately; otherwise, the 
link-layer will wait until new arrivals fill the queue with s min packets, at which instant a new 
transmission containing s min packets is scheduled. 

Let q n denote the queue occupancy at the n-th departure instant. Then q n evolves according 
to the following recursion: 

q n = min {V[s(g re _i)] + [<? n _i - s max ] + ,if - s(g„_i)} (7) 

where the notation [x] + represents the non-negative part of x, and V[s(q n ^ 1 )j is the number of 
packet arrivals during the transmission phase of the (n — l)-st inter-departure epoch, i.e., during 
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Fig. 4: Embedded Markov Chain at departure instants for the three possible different regions stated in ((6): dashed 
states are not reachable after a transmission starting at state i. 



the time required to transmit s(g„,-i) packets, from the moment they are scheduled until they 
are purged from the queue. Therefore, V[s(q n -i)} + [q n -i — s max ] + is the number of packets that 
would be present in the queue at the end of a departure if the buffer has infinite capacity, and 
K — is the maximum possible queue occupancy immediately after a departure. 

In order to find analytical formulation for key performance metrics such as delay and through- 
put, the steady-state queue occupancy probabilities need to be calculated. To derive these prob- 
abilities, first, the steady-state distribution for the queue states immediately after departures, 
7r d , is derived using the discrete-time embedded Markov chain method MIDI . Then the PASTA 



(Poisson Arrivals See Time Averages) property of the Poisson arrivals fllOl is applied to find the 
occupancy distribution at arbitrary times, 7r s , as a function of 7r d . 

In what follows, in order to make a clear distinction between the aforementioned two different 
steady-state probabilities, we define two different sets of states and a corresponding terminology 
for their probability distribution as follows: 

• queue state at departure instants: number of packets stored in the queue immediately after 
a departure. Hereafter, the steady-state probability distribution for these states is referred 
to as the departure distribution. This is what was denoted above by the row vector 7r d . 

• queue state at arbitrary times: number of packets stored in the queue at any arbitrary time. 
The steady-state probability distribution for these states will simply be referred to as the 
steady-state distribution of the queue. This is what was denoted above by the row vector 7r s . 

5.1.1 Departure Distribution, ir d 

The departure probability distribution, 7r d , is computed by solving the discrete-time embedded 
Markov chain§ of the occupancy of the batch-service queue immediately after departure instants, 
which can be done by solving the linear system 7r d = 7r d P, together with the normalization 

2. The Markov property of g„ is easily deduced from 0. 



11 



condition 7r d l T = 1, where P is the probability transition matrix of the embedded Markov 
chain, with each element p it j, i,j G [0, K — s m - m ], representing the probability to move from state 
i to state j. In this chain, transitions occur at the departure moments, and states represent the 
queue occupancy immediately after a departure. 

The irreducible embedded Markov chain consists of states i G [0,K — s m - m }. Note that, even 
though the queue has a capacity of K packets, a departure can never leave the queue in any state 
% > K — s min because at least s min packets are sent at every departure instant and these packets 
are held in the queue until their departure. Each transition occurs at a departure instant, which 
is immediately after receiving the last ACK, for a successful space-batch departure (Figure [5), or 
at the moment a packet is discarded, for the packets that reach the retransmission limit. Using a 
similar reasoning, from any given state i, the system can only transition into a subset of states, 
i.e., j G [[i - s max ] + , K - s(i)}. 

Figure H] depicts this Markov chain and the possible transitions from an initial state i for the 
three possible different cases depending on the value of v. 

1) When % < s min (Figure ©a)), the queue does not schedule a new transmission until there 
are at least s min packets available. Therefore, s(i) = s min for any i < s min , and transition to 
any state j is equivalent to transition from state s min to state j, and therefore, p it j = p Srnin ,j- 
In this case, the transitions are possible only into states j G [0, K — s min ]. 

2) When s min < i < s max (Figure HKb)), we have s(i) = i and therefore, transitions from state i 
can reach any state j G [0, K — i], 

3) Finally, when i > s max (Figure Htc)), we have s(i) = s max , and only transitions to states 
j G [i - s max , K - s max ] are possible. 

The probability of reaching a state j from any state i is given by: 



where j G [[i — s max ] + , K — s(i)]. For all other values of i, j, we have p iy j = 0. Note that j = K—s(i) 
is the last reachable state from any state i. Departing at state j = K — s(i) means that the queue 
has been containing q = K packets just before the departure, and therefore, some arrivals have 
possibly been blocked. For all other reachable states, the queue has had room for more packets 
just before the departure and therefore no arrivals could have been blocked. 

The number of arrivals during the transmission of a given space-batch depends on the time 
required for the transmission of the space-batch, which in turn depends on the space-batch 
size and the specific channel access protocol considered, as will be discussed in Section 15.21 
Let X(s(i)) be the random variable representing the service time for a space-batch of size s(i) 
packets, and fx(s(i))(x) be the probability density function (pdf) of X(s(i)). Note that X(s(i)) 
includes the time required for possible retransmissions due to collision or errors. 

For Poisson arrivals of rate A, the number of arrivals during X(s(i)) has the following distri- 



Pr {V[s(i)} = j — [i — s max ] + } , j < K - s(i) 
Pr {V[s(i)] > j — [i — s max ] + } , j = K - s(i) 



(8) 



bution: 




(9) 



12 



Ret. caused by a collision 

(packet arrival, the queue state q includes the new packet) 

»- q=1 q=2 q=3 blocked q=3 /' blocked 

COLLISION 



Ret. caused by tx. errors 



BUSY 



BO 

1 
DIFS 



ACKJout 
ACK * BO 



Tcf(1) 



MX 



X(s(1)) 



interdeparture epoch, T(1) 



q=0 



DIFS 




packet ACK tout 



blocked 

ACK tout 



AC K tout 



BO 



(errors) 



ACK 



BO packet i ACK 



Tc(2) 



" DIFS 



Tcf(2) 



DIFS 



Tcf(1) 



X(s(2)) 



interdeparture epoch, T(2) 



q=2 



(space-batch departure, the queue state q is just after the departure of the transmitted packets) 



Fig. 5: Temporal evolution of the queue state for a single node with a queue length of K = 3 packets, M > 2 
antennas, s m - m = 1 and s max = 2. After the departure of the first space-batch at state q = 0, the system remains 
in the empty state until a new packet arrives, at which instant a new space-batch transmission involving a single 
packet is initiated. The duration of this interdeparture epoch is T(0) = 4 +X(s(l)). The system departs at state q = 2 
since during X(l), three packets have arrived to the queue. Notice that the last arrival is blocked as it observes 
the queue fully occupied (q = K). A new space-batch is scheduled immediately at the beginning of the second 
interdeparture epoch, involving the two packets waiting for transmission. In this case the entire interdeparture 
epoch is spent in transmission mode and therefore, T(2) = X(s(2)). Observe how the second interdeparture epoch 
includes retransmissions caused by a collision (all packets are retransmitted) and transmission errors (only the 
non-ackonwledged packets are retransmitted). 



which in the particular case of exponentially distributed X(s(i)) with /ij = 1/ ' E[X(s{i))], it can 
be simplified to: 



Pi{V[s(i)] = v} 



Hi + A + A 



A 



(10) 



For any feasible state pair i.e., i G [0, A" — s min ] and j e \[i — s max ] + , K — s(i)], from © 



and (Q, we have: 



Pi,3 



(\f)v 

e- xt ^-f x{s{{)) (t)dt, j<K-s(i) 



o 



K-s(i)-\ 
z=i—s(i) 



(11) 



j = K - s(i) 



where v=j-[i- s max ] + . 

Finally from the departure distribution, the probability to transmit a space-batch of exactly s 
packets (at its initial attempt) can be computed by: 



r o 



$(s) = < 



E s min d 
9=0 71 9 

9=Smax 9 



S ^min 
S G ["S m i n ~\~ 1, S max 1] 
5 Smarx 



(12) 



5. 1.2 Steady-State Distribution, tt s 

Using the PASTA property IITOl of Poisson arrivals, the probability that at an arbitrary time in the 
steady-state the queue contains q = k packets is equal to the probability that a random arrival 
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observes k packets in the queue. In other words, 

irt = Pr{q a (t) = k} (13) 

where 7r| is the k-th element of 7r s , and q a (t) is the state of the queue observed by an arrival at 
time t. The right hand side of (fl~3~)) can be expanded by conditioning on q d {t), the state of the 
queue at the most recent departure before t, i.e., 



Til 



Pr {&(*) = k\q d (t) = 1} Pr {q d (t) = i} (14) 



The probability Pr{g d (£) = i}, can be viewed as the probability that the arrival at time t 
happens to occur during the departure state i of the embedded Markov chain discussed in 
the previous subsection. Therefore, this probability is equal to the expected fraction of time that 
the node spends in the departure state %, i.e., 

p,fa W =o= ; ? T (i)1 as) 

where irf is the j-th element of n d , the departure distribution, T{%) is the random variable 
denoting the time spent in departure state i, and therefore E[T] is the expected length of an 
interdeparture epoch. 

To calculate E[T(i)] note that of the time spent in departure state i, X(s(i)) seconds will be 
spent in transmission mode, and if i < s min , an additional I(s(i)) seconds will be spent in idle 
mode before entering transmission mode. Therefore: 

E[T(i)\ = E[I(8(i))]+E[X(8(i))] 

= j[s min -i} + + E[X(s(i))} (16) 

where j- [s min — i} + is nonzero only when % < s min and is equal to the expected time needed 
for the queue occupancy to reach s min . An example of the temporal evolution of the system is 
depicted in Figure |5j 



The term Pr {q a (t) = k\qd(t) = %} in (|T4)) is the probability that an arrival during the departure 
state i observes k packets in the queue, which can be viewed as the expected fraction of arrivals 
in departure state i that observe k packets in the queue. The expected total number of arrivals 
in state i is given by XE[T(i)]. Of these, only one may observe k < K frames in the queue, 
provided that there are enough arrivals. Let q n+1 be the state at which the next departure will 
leave the queue. Then the queue occupancy just before this next departure is q n+l +s(i). In order 
for an arrival to have observed k packets in the queue, we need q n+ \ +s(i) > k + 1. The expected 
number of arrivals in state i that observe k packets in the queue is given by: 

K-s(i) 

Pr {q n +i >k + l- s(i)\q n = i} = ^ Vi,j (17) 

j=k+X— s(i) 
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Based on this, we have 

Pr {&(*) = k\q d (t) = *} = Ej 7t^; ( f tJ (18) 

XE[T{i)\ 



From (dH), ((15|) , and (|T8l> , the steady state queue occupancy distribution, 7r s , for states < k < 



K — 1 can be computed as 

A^E-q E , o< fc <A--i 

vr| = <( L x _;=° \i=*+i-w / (19) 

i - *?> & = * 

i=0 

5.2 Medium Access Protocol 

As we saw in the previous subsection, in order to find the steady-state queue occupancy dis- 
tribution, we need to calculate E[X(s)}, the expected service time for a space-batch of size s 
packets. Since all the nodes are assumed to have the same behavior, when calculating different 
MAC layer parameters, we will use this symmetry and focus on a single reference node. 

The time required for the successful transmission of a space-batch depends on the number 
of times it has to be retransmitted due to collisions or transmission errors, the duration of 
each transmission, which itself depends on the number of packets sent, and finally the random 
duration of the backoff pauses caused by transmissions from other nodes. Therefore, in order 
to calculate E[X(s)], we need the probability that the reference node collides with any other 
node, the packet error probability, and the probability that any other node transmits while the 
reference node is in backoff, and the average duration of the pauses caused by the latter. The 
probability that a transmitted packet contains errors was previously discussed and is given by 
fl]}. In this section, we turn our attention to the Medium Access Control to calculate the rest of 
the aforementioned probabilities. 

We adapt the definition of time slots used in llT8i which defines a slot to be the time in 
between two decrements of the backoff counter, which is not constant, as the backoff counter is 
paused every time a transmission is detected on the channel. 

5.2. 1 Conditional Collision Probability 

Under saturated conditions, a node will always have packets to transmit, and therefore, after 
completing every transmission, it will select another random backoff value and start counting 
down. In this case, the node's backoff counter reaches zero on average every E[B] + 1 slots, 
where E[B] = cw 2 ~ 1 is the expected number of slots a random backoff will last. Therefore, the 
probability that a node in saturation transmits in a randomly selected slot in steady-state is 
given by r s = 

When the network is not saturated, a node can transmit only if it has at least s min packets 
ready to be transmitted and its backoff counter has reached at the end of the previous slot. 
Therefore, for unsaturated traffic, the probability that during a given backoff slot of the reference 
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node, another given node starts transmitting can be approximated by [19J: 

T = pT * = Efiri (20) 

where p is the steady-state probability that the node has at least s min packets stored in the queue, 
p — 1 — Xli=o n_1 n i- Note that under saturated conditions, p = 1 and therefore, (|20]> can also be 
used for saturated traffic. 

Since the nodes sense the channel at the beginning of every backoff slot and pause their 
backoff if they sense a transmission, a transmission from the reference node will collide only if 
at least one other node starts transmitting at the same time as the reference node. Therefore, the 
conditional collision probability, i.e., the probability that a transmission by the reference node 
ends up in a collision is given by: 

p = 1 - (1 - r)^- 1 (21) 
5.2.2 Expected Number of Retransmissions 

Space-batch retransmissions are caused by collisions or transmission errors. In case of a collision, 
all packets included in the space-batch are lost and should be retransmitted. In case of a collision- 
free transmission containing m packets, from to m frames can contain errors, each one with 
probability PER(m), independently of the rest. In this case, only the unacknowledged packets 
will be retransmitted. Let T c f(s) be the average number of collision-free transmission attempts 
required to successfully transmit all the s packets included in the initial space-batch, which in- 
cludes all the erroneous transmission attempts and a single successful one. Let T c be the average 
number of transmission attempts required for every collision-free transmission, i.e., the number 
of collisions before every collision-free transmission plus the collision-free transmission itself. For 
example, CE|E|CCE|CS is a possible transmission attempt sequence outcome for transmitting a 
given space-batch, with C, E, and S indicating collided, erroneous, and successful transmission, 
respectively. In this example, the number of collision-free (but possibly erroneous) transmission 
attempts is 4. The last collision-free transmission attempt is the only one that is successful. Then, 
the total number of transmission attempts required per space-batch is T(s) = T c f(s)T c . 

Assuming that the collision probability, p, remains constant regardless of the number of 
collisions that the target space-batch has already suffered, the number of attempts per collision- 
free transmission follows a geometric distribution, and its average value, T c is given by: 

Tc-V^l-p)/-^- (22) 

To compute T c f(s), we need to keep track of the unacknowledged packets remaining for 
retransmission, for which we use an absorbing, time-homogeneous Markov chain, whose state 
space is given by S = {m}^ =0 . Each state m represents the remaining number of unacknowledged 
packets, and s, the original size of the space-batch. The Markov chain reaches the absorbing state 
m = when all the packets have been successfully transmitted. Transition from a state i to a 
state j requires the successful transmission of i — j of the i packets transmitted in state i, which 
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happens with probability: 



V 



i,3 



i-j 



[1 - PER(z)) i_i PER(z) j 



(23) 



Let P be the transition probability matrix of the chain, with at its element. Given 
that the state (0, 0) is absorbing, the first row of P contains all zeros except for its first element 
which is one. In this case, P can be represented in the following canonical form: 



TR 



ABS 



/ 



TR 


ABS 




O(lxs) 


R-(sxl) 


Q(sxs) 



(24) 
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where the rows and columns corresponding to the transient and absorbing states are indicated 
by TR and ABS, respectively. For an absorbing chain with the canonical representation above, 
the following properties hold 1120 II : 

1) The matrix I — Q has an inverse given byN = I + Q + Q 2 H — • , with each of its elements, 
riij, being the expected number of times the chain visits state j before absorption, given 
that it started in state i. 

2) Let T C f(z) be the expected number of transitions before absorption, starting from state i, 
and T c f a column vector whose z-th entry is T c f(z). Then Y c f = N ■ c, where c is an s x 1 
column vector containing all ones. 

The number of required transmissions to successfully transmit all packets in a space-batch 
initially containing s packets, T C f(s), is equal to the expected number of steps to absorption in 
the aforementioned absorbing Markov chain, starting at state s. By solving the chain only once 
for s = M, we can obtain all T c f(i) values for all i. 



5.2.3 Expected Service Time 

The expected service time for a space-batch of size s, E[X(s)], is the time since it is scheduled 
until either the successful arrival of all of its packets is confirmed, or it is discarded due to 
reaching the maximum number of retransmissions, and is given by: 

E[X(s)] = T(s)E[B] 1 + T rf (s) (T d (s) + (T c - l)T c (s)) (25) 

where 7 is the average duration of a backoff slot and includes possible backoff countdown 
interruptions, and T c f(s) and T c (s) are the expected duration of collision-free and collided 
transmission attempts, respectively, during the transmission of a space-batch of initial size s. 
The first term in the right-hand-side of ((251) , represent the total amount of time spent in backoff 
over all transmission attempts required to send a space-batch of size s, and the second term the 
time spent in transmission and receiving acknowledgements. The expected service time averaged 
over the space-batch size, s, is simply E[X] = Es=sl ln &( S ) E [ X ( S )}- 
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As mentioned before, the duration of a slot can vary depending on whether or not at the 
beginning of that slot a transmission is detected on the channel. Let a be the duration of an 
empty slot, i.e., the time between two backoff decrements if the channel is detected idle. Then 
the average slot duration, 7, can be calculated as: 

7 = Pe a + p cf (r rf + DIFS + a) + Pc (T c + DIFS + a) (26) 

where p e , p c {, and p c are the probabilities that during a given backoff slot of the reference node 
the channel is empty, contains a collision-free transmission, or contains a collision, respectively. 
T c f and T c are the average duration of a collision-free and a collided transmission observed by 
the reference node in backoff. The DIFS and extra slot duration added to T C f and T c account for 
the fact that in both cases of collision-free transmission and collision, the reference node has to 
sense the channel idle for at least a seconds before decrementing its backoff counter and that 
the channel is only detected idle after sensing it free for DIFS seconds. 

A node performing backoff will find the channel idle if there is no other node transmitting 
at that slot. Therefore, the probability p e , that a given node finds the channel empty during a 
backoff slot is given by: 

p e =(l-r) iV - 1 (27) 

The probability of finding a successful transmission during a backoff slot, pa, is equal to the 
probability that only one of the other N — 1 active nodes transmit during that slot, i.e., 

p c{ = (N-l)-T(l-rf- 2 (28) 

Finally, the probability that a selected slot contains a collision, p c is the complementary of the 
two previous cases 

Pc = 1 - Pcf - Pe (29) 

To compute T cf and T c it should be noted that both are expected durations averaged over all 
possible values of m, the number of packets included in each transmission. The number of pack- 
ets included in a transmission attempt depends on the initial number of packets, s, scheduled for 
that space-batch and the number of packets not successfully transmitted in previous attempts. 
Let p m \ s denote the probability that a given transmission attempt during the transmission of a 
space-batch of initial size s, happens to contain m packets. Then, the average duration for each 
transmission attempt for a space-batch that of initially s packets, T c f(s), is given by: 

s 

T c f(s) = ^Pm\ s Td(m) (30) 

m=l 

where T cf (m) is given by ©. Averaging for different values of s results: 

T cf = *(s)Tcf(s) (31) 

S — S m i n 
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TABLE 1: Specific parameters considered for the MACMA space-batch and ACK frames 



Field 


Length 


Field 


Length 


Preamble 


128 bits 


SFD 


16 bits 


SIG 


48 bits 


Training Sequence 


64 bits 


FC 


16 bits 


Tx Address 


48 bits 


SC 


16 bits 


Rx Address 


48 bits 


Packet length 


variable 


CRC 


32 bits 



Using a similar approach, for T c , we first calculate the duration of a collision given that one 
of the two nodes involved in the collision has transmitted a space-batch of initially s frames: 

Smax Si S 

Sl=S m j n 1=1 j = l 

where to simplify, we have assumed that the probability of having more than two nodes involved 
in a single collision is negligible. Averaging over all possible s values, we get: 

T c = *(s)T c (s) (33) 

S — Smin 

The conditional probability of transmitting a space-batch containing m packets when the initial 
space-batch included s packets, p m \ s , is calculated using a Markov chain with state space S = 
{mY m=1 , where s is the size of the initial space-batch. This Markov chain is exactly the same as 
the absorbing Markov chain used to calculate T c f(s), except that, in order to be able to compute 
the stationary distribution, the absorbing state has been removed. In this case, when all packets 
have been successfully transmitted, the chain goes back to the initial state s, accounting for the 
next scheduled space-batch of initial size s. 

Let pij be the probability to move from a state i to a state j, given by: 

{Pij 1 < i < s, 1 < j < i 

Pi,o 1 < i < s, j = s (34) 

Ps,0 + Ps,s i = s, j = s 

where pij are given by (|23~1> . 

Let P s be the transition matrix of this Markov chain and tt s , the vector containing its stationary 
probabilities. Then, tt s is obtained by solving the equation system tt s = tz s P s , together with the 
normalization condition 7r s l T = 1. Then, p rn \ s is the m-th element of iz s , i.e., p m \ s = 7r s (m). 

6 Performance Evaluation 

In this section, the performance of the MACMA protocol is evaluated through simulation and 
compared to the results obtained from the analytical formulation presented in the previous 
sections. 
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6.1 Scenario, Simulation, and Model Parameters and Considerations 

In order to evaluate the performance of the MACMA protocol and validate the presented 
analytical model, we consider the scenario described in Section [3j A simulator of this scenario has 
been built in C++, from scratch, based on the COST (Component Oriented Simulation Toolkit) 
libraries El . The data transmission rate and the physical transmission rate are set to r = 100 and 
?>hy = 10 Mbps respectively, and a constant packet length of L = 8000 bits are used. Additionally, 
for the MAC parameters, a single-stage backoff with CW = 32, DIFS = 50 fj,s, SIFS =10 fj,s, and 
er = 20 /is is considered. Each node has a queue of size K = 50 packets. In Table [TJ the values for 
different fields of MACMA space-batch and ACK frames are shown (Figure [3]). All the values 
considered are based on those used in |;3j. With respect to the PER, we have considered a 
reference SNR of £ re f = 15 dBs. 

Regarding the analytical model, when numerically evaluating the formulation derived in the 
previous sections, the service time is assumed to follow an exponential distribution, fx(s(i)) — 
Hie'^ 1 , where /ij = l/E[X(s(i))]. This is done as a tradeoff between accurately capturing the vari- 
ations in the service time caused by the channel access mechanism, and reducing the numerical 
complications introduced when solving the model using more accurate distributions (which are 
different for each case and not known a priori). Moreover, even though it has been proven that 
the service time for the DCF in both saturated and non-saturated conditions follows a skewed 
distribution ||22|, the assumption of exponentially distributed service times has been widely used 
as it provides a reasonable level of accuracy ||22|. However, if a more accurate model is required, 
the specific service time distribution can be computed and included in the model in a similar 
way as described in 1123 II , llTil . 

Two performance metrics are considered for the evaluation of the MACMA protocol: (1) ag- 
gregate throughput (bits/second) is defined as the number of useful data bits that are successfully 
transmitted in the network per unit of time and is given by S = N(l — p\,)\L, where pb = Tffc 
is the blocking probability, i.e., the probability that a packet is lost due to finding a full queue 
upon its arrival, and (2) average queueing delay (seconds) is defined as the amount of time a 
packet spends waiting in the queue, from the moment it is admitted to the queue until it is 
purged from the queue after being successfully transmitted or discarded. Using Little's law, the 
average queueing delay is given by E[D\ = X (i-l h ) > wnere E[Q\ = Y1^=q Q^l' is the average queue 
occupancy. 



6.2 Results 



Figures 6(a) and 6(b) show the aggregate throughput and average delay of MACMA for a 
different number of nodes when £ = 20 dB, s min = 1, s max = 2, and M = 2 or 4 antennas. 
With respect to the aggregate throughput (Figure |6(a)[ ) and the average delay (Figure |6(b)| >, the 
gain observed when the nodes are equipped with M = 4 antennas is due to the low PER that 
the Zero Forcing receiver is able to provide by using the extra channel's degrees of freedom for 
diversity [9j. A similar gain is observed in both the aggregate throughput (Figure |6(c)| > and the 
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Fig. 6: The impact of s max on the performance of MACMA in non-saturation conditions with M = 4 antennas. In 
the legend, each curve is marked by the value of [M, s m ; n , s max ]. The x-axis and y-axis label 'SNR' in the figures 
refers to £o- 



average delay (Figure 6(d)| > when the number of nodes is fixed to iV = 6 and 12 nodes, and the 
packet arrival rate at each node increases from 62.5 to 1000 packets /second. 

The impact of different s m i n values is evaluated when multiple nodes contend for the channel. 
The results are plotted in Figure [7J The considered parameters are M = 4 antennas, £o = 40 dBs, 
Smax = 4, and s min = 1 or 4 packets. Setting s min = 4 extends the performance of the MACMA 
protocol, supporting a large number of nodes compared to the case with s mm = 1 (Figure |7(a)| >. 
Additionally, an interesting issue here, already known from the analysis of the DCF with single- 
antenna nodes [24], is that in some specific circumstances (i.e., when the collision probability 
just before saturation remains very low), the non-saturated throughput just before the saturation 
point is higher than the saturation throughput, as it is also observed here. Setting s min = 4 results 
in fewer channel attempts as the queue remains inactive, i.e., with less than s min packets, for 
longer periods of time (Figure 7(b)| >, which results in a lower collision probability until the 
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Fig. 7: The impact of s m ; n on the performance of MACMA in non-saturation conditions. In the legend, each curve 
is marked by the value of [M, s m i n , s max ]. The y-axis label 'SNR' in the figures refers to £o- The traffic load in each 
node is 2 or 4 Mbps. 



saturation point is reached (Figure |7(c)[ ). However, values of s min larger than 1 have a negative 
impact on the delay (Figure 17(d)) ), due to the extra time required to have s min packets in the 
queue at low traffic loads (small N). Note that this time is proportional to J2i=f j( s min — i) 
seconds. 



7 Conclusions 

In this paper, an analytical model for a CSMA/CA (DCF-like) MAC protocol with MPTR capa- 
bilities, called MACMA, was presented. The model was built by combining a batch-service 
queueing model with a non-saturated model of the DCF Simulation results show that the 
presented model is able to capture the fundamental interactions between the CSMA/CA channel 
access mechanism and the MPTR capabilities of the MACMA protocol, as well as providing a 
very accurate estimation of the system performance. 
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As it was observed, the required extra overheads to support multiple packet transmission 
are compensated by the larger amount of data that can be sent at each transmission, as well 
as by the reduction in the number of transmissions, which in turn results in lower collision 
probabilities for a given number of contending nodes and a given traffic load, thus providing a 
higher throughput and lower delay 

Moreover, given a node equipped with multiple antennas, two new parameters called s min 
and s max , which control the number of packets that are scheduled at each transmission, were 
introduced and their impact on the system performance was evaluated, drawing several con- 
clusions. Mainly, it was observed that, although the number of packets that can be included in 
a space-batch increases when s max increases, the PER also increases. Therefore, the value of s max 
has to be set to the highest possible value that still guarantees a reasonably low PER (i.e., based 
on the observed SNR and available number of antennas in each node). It was also observed that 
adapting s min to the traffic load and the number of active nodes, the system performance can 
be optimized in terms of delay and throughput. Moreover, although not analyzed in this paper, 
the impact of both parameters also depends on the backoff algorithm and its parameters (i.e., 
CW and number of stages), thus opening a new opportunity to search optimal configurations 
for all those parameters simultaneously. 
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